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1. Introduction 

At present, most lattice QCD simulations are performed at unphysical light dynamical quark 
masses. Fitting of lattice data to forms calculated in chiral perturbation theory (XPT) [jlj, ^ makes 
possible a controlled extrapolation of lattice results to the physical light quark masses and to the 
chiral limit. This approach also allows one to determine the values of low-energy constants (LECs) 
in the theory, which are of phenomenological significance. Although three-flavor XPT has been 
used successfully for simulations with 2+1 dynamical quarks, we are still interested in the applica- 
tions of two-flavor XPT for the following reasons: 

1. The up and down dynamical quark masses in simulations are usually much smaller than the 
strange quark mass, which is near its physical value, hence SU(2) XPT may serve as a better 
approximation and probably converges faster than SU(3) XPT. 

2. Fits to SU(2) XPT can give us direct information about the LECs in the two-flavor theory, 
especially I3 and I4. 

3. By comparing results from these two different fits, we can study the systematic errors result- 
ing from the truncations of each version of XPT. 

Recently, some groups have used SU(2) XPT for chiral fits to data from three-flavor simulations [||, 
IJ. Here, we perform such an SU(2) chiral analysis for MILC data from simulations with 2+1 
flavors of staggered fermions. 



2. Rooted SU(2) SZPT 

For staggered quarks, the correct effective field theory is staggered chiral perturbation theory 
(SZPT) []5| |[ 0, H ||], in which taste- violating effects at finite lattice spacing are incorporated 
systematically. Physical quantities expressed in SXPT become joint expansions in both the quark 
mass m q and a 2 , where a is the lattice spacing. 

For each quark flavor, there are four species (tastes) in the continuum limit. To obtain physical 
results, we use the fourth root procedure to get a single taste per flavor in the continuum limit. 
Although it has been shown that this procedure produces violations of locality at non-zero lattice 
spacing non-perturbatively [|K]], recent work indicates that locality and universality are restored in 



the continuum limit. For a recent review of the fourth-root procedure see Ref. flllQ and references 
therein. 

In the two-flavor case, only up and down quarks appear in the chiral theory. Correspondingly, 
there are only pions, and no kaons, in SU(2) SZPT The staggered Lagrangian is formulated in 
the same way as in Ref. [g], except that those parts related to the strange quark are omitted. Fol- 
lowing the procedures used in the three-flavor case, one can calculate the partially-quenched light 
pseudoscalar mass and decay constant through NLO. The result is [|I"2|]: 
2 
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where I3 and I4 are the standard SU(2) XVY LECs, and p\ and P2 are two extra NLO LECs that 
enter in the partially-quenched case. 8y, 8' A are taste- violating hairpin parameters, and L'^, L'L^ 
are taste-violating analytic LECs. Chiral logarithms l(m 2 ), l(m 2 ) and residue functions R, D are 
given in Ref. [H], with the denominator mass-set arguments in the SU(2) case defined as: 



[21 

xy,} = {m Xn m Yl }, 



i-^XYy } = i m Xv . m Yv 



(2.3) 



The numerator mass-set arguments of the residues are always {/is} = {wiy.}, where the taste label 
S is taken equal to the taste of the denominator set. 

To Eqs. (2T) and ( |2.2| ), we add the NNLO chiral logarithms that were calculated by Bijnens 



and Lahde [|13j. Since taste splittings are not included at NNLO, there is an ambiguity in defining 
the pion mass in the continuum formulae. In practice, we use the root mean square (RMS) average 
pion mass in calculations of NNLO chiral logarithms. This is systematic at NNLO as long as the 
taste splittings between different pions are significantly less than the pion masses themselves. This 
condition is best satisfied on the superfine and ultrafine lattices. 



3. Ensembles and Data Sets 



At the present stage, we have the MILC data for the light pseudoscalar mass and decay constant 
at five lattice spacings from 0.15fm to 0.045 fm, generated with 2+1 flavors of asqtad improved 
staggered quarks. For each lattice spacing, we have several different sea quark masses as well as 
many different combinations of valence quark masses. In order for the SU(2) formulae to apply, 
we require both sea and valence quark masses to be significantly smaller than the strange quark 
mass, i.e., m s * a <^niK, and m™ lence ^vyik- In the fits described below, we use the following cutoff 
on our data sets: 

mi < 02m p s hys , m x + m y < 0.5m p s hys , (3.1) 

where m; is the light sea quark mass, and m x and m y are the valence masses in the pion. 

To be able to consider the strange quark as "heavy" and eliminate it from the chiral theory, it 
is also necessary that taste splittings between different pion states be much smaller than the kaon 
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Ensemble 


ami 


am s 


P 


size 


m n L 


« 0.09 fm (F) 


0.0062 


0.031 


7.09 


28 3 x 96 


4.14 


« 0.09 fm (F) 


0.00465 


0.031 


7.085 


32 3 x 96 


4.10 


« 0.09 fm (F) 


0.0031 


0.031 


7.08 


40 3 x 96 


4.22 


» 0.09 fm (F) 


0.00155 


0.031 


7.075 


64 3 x 96 


4.80 


« 0.06 fm (SF) 


0.0036 


0.018 


7.47 


48 3 x 144 


4.50 


« 0.06 fm (SF) 


0.0025 


0.018 


7.465 


56 3 x 144 


4.38 


« 0.06 fm (SF) 


0.0018 


0.018 


7.46 


64 3 x 144 


4.27 


« 0.045 fm (UF) 


0.0028 


0.014 


7.81 


64 3 x 192 


4.56 



Table 1: Ensembles used in this analysis. The quantities ami and am s are the light and strange sea quark 
masses in lattice units; m K L is the (sea) Goldstone pion mass times the linear spatial size. The fine ensembles 
are not used in our central value fit, but only in estimating systematic errors. 



a 


» 0.09 fin (F) 


« 0.06 fm (SF) 


ps 0.045 fm (UF) 


ams 


0.031 


0.018 


0.014 


ami 


0.00155 


0.0062 


0.0018 


0.0036 


0.0028 


m K (MeV) 


574 


613 


525 


548 


565 


m Goldstone (MeV) 


177 


355 


224 


317 


324 


m™ s (MeV) 


281 


416 


258 


341 


334 


m^MeV) 


346 


463 


280 


359 


341 



Table 2: Kaon masses and lightest (sea) pion masses on some sample ensembles. Here three different pion 
masses are shown: Goldstone, RMS and singlet. r\ = 0.31 17 fm is used. 

mass. Furthermore, taste splittings should be significantly smaller than the pion mass itself for the 
continuum formulae for the NNLO chiral logarithms to be approximately applicable. 

The lattices that are at least close to satisfying all these conditions include four fine (a « 0.09 fm) 
ensembles, three superfine (a sa 0.06 fm) ensembles and one ultrafine ensemble (a w 0.045 fm). Rel- 
evant parameters for these ensembles are listed in Table |]. 

In Table ||, we list the Goldstone, RMS and singlet pion masses on representative ensembles. 
It can be seen that for the fine (a rj 0.09 fm) ensembles, either some pion masses are close to the 
kaon mass, as on ensemble (ami,am s ) = (0.0062,0.031), or the taste splittings between pions are 
comparable to the pion mass, as on ensemble (ami,am s ) = (0.00155,0.031). As a result, the data 
from fine lattices may not be well described by SU(2) formulae with continuum NNLO chiral 
logarithms. Our central fit uses superfine and ultrafine data only, while we include fits to all three 
kinds of lattices to estimate systematic errors. 

There are a total of 29 parameters in our fits. The following list shows how these parameters 
are treated in the central fit. 

(a) LO: 2 unconstrained parameters, pL and /. 

(b) NLO (physical): 4 parameters, I3, I4 and two extra LECs p\,p2 that only appear in 
partially-quenched XVT. All of these parameters are unconstrained. 
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(c) NLO (taste- violating): 4 parameters. 8y,8^ are constrained within errors at the values 
determined from SU(3) SZPT fits [14, 15]; L'L\ and L'^\ are constrained around 0, with width 
of 0.3 as estimated in Ref. @|. 

(d) NNLO (physical, 0(p 4 )): 5 parameters (h, 1%, h, p$, that first appear in meson 
masses and decay constants in the NNLO chiral logarithms. l\ and l 2 are constrained by the 
range determined from continuum phenomenology [16]; lj is not constrained since it is not 
directly known from phenomenology [jl^]. The partially-quenched parameters p$ and p4 are 
not constrained. 

(e) NNLO (physical, <^(p 6 )): 8 parameters Cf, constrained around with width 1 in "natural 
units" (see Ref. Q). 

(f) The physical LO and NLO parameters are allowed to vary with lattice spacing by an 
amount proportional to a s (aA) 2 , which is the size of the "generic" discretization errors with 
asqtad quarks, where A is some typical hadronic scale. This introduces 6 additional parame- 
ters that are constrained around with width corresponding to a scale A = 0.7 GeV. 

Alternative versions of the fits, in which the width of the constraints are changed, or some con- 
strained parameters are left unconstrained (or vice versa), have also been tried, and the results from 
those fits are included in the systematic error estimates. 



4. Preliminary Results 

For the central fit, we use three superfine ensembles {ami,am s ) = {(0.0018,0.018), (0.0025, 
0.018), (0.0036,0.018)} and one ultrafine ensemble (am u am s ) = (0.0028,0.014). This fit has a 
% 2 of 37 with 33 degrees of freedom, giving a confidence level CL = 0.3. The volume dependence 
at NLO has been included in the fit formulae. A very small (< 0.3%) correction for "residual" finite 
volume effects [|n|, [l8|] is applied at the end of the calculation and incorporated in the systematic 
errors of our final results. 

In Fig. |l| we show the fit results for f n and m 2 n / [m x + m y ) as functions of the sum of the quark 
masses (m x + m y ). The red curves show the complete results through NNLO for full QCD in the 
continuum, where we have set taste splitting and taste-violating parameters to zero, extrapolated 
physical parameters as a — ► linearly in a s a 2 , and set valence quark masses and light sea quark 
masses equal. Continuum results through NLO and at tree level are shown by blue and magenta 
curves, respectively. It can be seen that the convergence of SU(2) ZPT is much better for the decay 
constant than for the mass. Nevertheless, the chiral corrections in both cases appear to be under 
control. 

The SU(2) plots presented previously in Ref. jl9| ] are somewhat different from those shown 
here because the earlier fits allowed for a 2 variations in the NNLO analytic parameters (c ; ). Such 
variations are of higher order than the NNLO terms included in this work. 

At the last step, we find the physical values of the average u,d quark mass m by requiring that 
the 71 has its physical mass, and then find the decay constant corresponding to this point in Fig. |l] 
(left). With the scale parameter r\ = 0.318(7) fm [ p"4| ] determined from T-splittings, we obtain the 
result for f K : 

f % = 128.3(9) MeV (4.1) 
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f„ ; SU(2) xPT 



O.lf 



0.17 



[w 

> 0. 



16 



0.15 



0.14 



+super fine, anig=0.01B 




oultra fine, am^O.014 




- CL = 0.3 




X Z /d°f=37/33 




- yy^Y NNLO, cont. 


0.0036 - 
+ 0.0025 - 


///y — NLO, cont. 


+ 0.0018 - 
0.0028 - 


- w 

w * extrap 




- / ti expt. (r,=0.318fm from T) 




/ LO, cont. 


m K +m y =0.5my 


i i i i i i i i 1 i i 





0.00 0.02 

(m x +m y )r, x ,Z„ /.^ ) 

(a) 



0.04 0.06 

finc\ 

m/ m 



6.6 — 



: e 6.4 



£ 6-2 



6.0 



m=/(m x + m y ); SU(2) X PT 

♦ super fine, am;=0.018 CL = 3 
o ultra fine, am;=0.014 xVdof=37/33 

0.0036 
+ 0.0035 
+ 0.0018 
O 0.0028 _ 

. LO, cont. 




0.00 0.02 

K+m ; )r, x (7 m 7,r) 

(b) 



0.04 0.06 

fine"\ 
m/ '-'m 



Figure 1: SU(2) chiral fits to /^(left) and m|/ (m x + m,,)(right). Only points with the valence quark masses 
equal (m x = m y ) are shown on the plots 



where the first error is statistical and the second is systematic. This agrees with the PDG 2008 
value, f n = 130.4 ± 0.2 MeV Alternatively, using the pion decay constant from NNLO SU(3) 

0.3117(6) 



XPT fits to define our scale gives r\ 

fz = 123.7(9) (18) MeV 
h =3.0(6) 
m = 3.21 (3) (5) (16) MeV 



-31 ; f m Pl^O- With this new r\, we obtain: 

B2 = 2.89(2) (lg) (14) MeV 
f 4 = 3.9(2)(3) 

(uu)z = -[280(2) (4) MeV] 3 



(4.2) 



The quark masses and chiral condensate are evaluated in the MS scheme at 2GeV. We use the 



two-loop renormalization factor in the conversion Q21p. Errors from perturbative calculations are 
listed as the third error in these quantities. All the quantities agree with results from SU(3) SZPT 



fits [15] within errors. 



5. Discussion and Outlook 

We have performed NNLO SU(2) chiral fits to recent asqtad data in the light pseudoscalar 
sector. Results for SU(2) LECs, the pion decay constant, and the chiral condensate in the two-flavor 
chiral limit are in good agreement with those obtained from NNLO SU(3) fits (supplemented by 
higher-order analytic terms for quantities involving strange valence quarks) Jl5|]. It can be seen 
from our plots that SU(2) XVI within its applicable region converges much faster than SU(3) XVI. 
For the point 0.05 on the x-axis in Fig. [I], the ratio of the NNLO correction to the result through 
NLO is 0.3% for f n and 2.6% for m K /(m x + m y ). In contrast, the corresponding numbers in the 



SU(3) fits are 2.9% and 15.6% respectively (Fig. 2 of Ref. [15]), although the large correction in 
the mass case is partly the result of an anomalously small NLO term. Note that the SU(3) plots 
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use a non-physical strange quark mass, m s = 0.6m% ys , while for the SU(2) plots, the strange quark 
mass is near the physical value, m s pa m p s h} J . This explains why the two-flavor chiral limits on the 
SU(3) and SU(2) plots are not the same. 

Since the simulated strange quark masses vary slightly between different ensembles, the pa- 



rameters in SU(2) SZPT should also change with ensemble Q12[]. We plan to incorporate this effect 
in our fit to see if we can improve the confidence levels. Another step would be to include the kaon 
as a heavy particle in SU(2) SZPT [0] in order to study the physics involving the strange quark, 
e.g., the kaon mass and decay constant. This approach has recently been used in Refs. 
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